An Infrared Safe perturbative approach to Yang-Mills correlators 
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We investigate the 2-point correlation functions of Yang-Mills theory in the Landau gauge by 
means of a massive extension of the Faddeev-Popov action. This model is based on some phe- 
nomenological arguments and constraints on the ultraviolet behavior of the theory. We show that 
the running coupling constant remains finite at all energy scales (no Landau pole) for d > 2 and 
argue that the relevant parameter of perturbation theory is significantly smaller than 1 at all ener- 
gies. Perturbative results at low orders are therefore expected to be satisfactory and we indeed find 
a very good agreement between 1-loop correlation functions and the lattice simulations, in 3 and 4 
dimensions. Dimension 2 is shown to play the role of an upper critical dimension, which explains 
why the lattice predictions are qualitatively different from those in higher dimensions. 
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I. INTRODUCTION 

There is up to day no fully satisfying covariant gauge 
fixing of non-abelian gauge theories. The common 
Faddeev-Popov (FP) procedure is known to be invalid be- 
cause of the Gribov ambiguity [ij: the gauge constraint 
(such as = in the Landau gauge) have many 

solutions (Gribov copies) equivalent up to a gauge trans- 
formation. Taking into account this ambiguity would 
give non-perturbative contributions (typically of the form 
exp(— cst./g^)) that vanish at all orders of the perturba- 
tion theory. Consequently, the FP procedure is actually 
justified when studying high energy phenomena and the 
perturbative predictions are in excellent agreement with 
the experiments. On the other range of the spectrum, 
when considering infrared (IR) properties, the perturba- 
tive analysis extracted from the FP procedure is incon- 
sistent because the coupling becomes large. The coupling 
is found to diverge at a finite energy scale (known as a 
Landau pole) but of course this prediction is out of the 
range of validity of perturbation theory. There are two 
possible explanations of the failure of the perturbative 
predictions of the FP procedure. Either the coupling in- 
deed reaches large values (but this seems to be at odds 
with the lattice simulation results) or some nonpertur- 
bative effects, such as the Gribov ambiguity, invalidates 
the FP procedure itself. 

In fact, we can avoid fixing the gauge by simulating 
gauge invariant quantities on the lattice and most of the 
IR studies are done in this way. However it would be 
convenient to have some analytic (or semi-analytic) pre- 
dictions in that regime, which seem to require to fix the 
gauge (see however 0-i4!|), and therefore to take into ac- 
count the Gribov ambiguity. 

Several analytic methods have been considered in or- 
der to access the IR properties. The most developed 
is based on the Schwinger-Dyson (SD) equations E5:-16[, 
which consist in an infinite set of coupled equations for 
the vertex functions. In order to make predictions, it 



is necessary to truncate in some way this infinite set. 
Different schemes have been proposed but most of them 
consider the equations for the 2-point functions with an 
ansatz for the 3 and 4-point functions. (Most of the anal- 
yses have been done in the Landau gauge and we restrict 
the discussion to this gauge in the rest of the paper.) This 
leads to some predictions on the behavior of the ghost and 
gluon propagators and two types of solutions have been 
found. In the so-called scaling solution 

both 

propagators have a power-law behavior in the IR, the 
ghost propagator is more singular than the bare one and 
the gluon propagator approaches zero in that limit. The 
exponents governing these two power-laws are not inde- 
pendent. In the so called decoupling solution (or massive 
solution) the gluon propagator goes to a con- 

stant in the IR and the ghost propagator is as singular 
as in the bare theory. 

These correlation functions were also studied in the 
framework of the Non-perturbative Renormalization 
Group (NPRG) equation Again, one has to make 

some truncations in an infinite set of coupled equations 
[H, [S 1131 • The results again depend on the approxima- 
tion scheme, but seem to be consistent with the scaling 
solution. 

A third approach, known as the Gribov- Zwanziger 
(GZ) model Q, [SI, [H, relies more specifically on the 
influence of the Gribov copies. By introducing several 
auxiliary fields, it is possible to restrict effectively the 
functional integration on the gluon field to the first Gri- 
bov region where the FP operator is positive definite. It 
was originally expected that this restriction would give 
an unambiguous gauge fixing but it was later realized 
that many copies of some field configuration were present 
within the first Gribov region This procedure sig- 

nificantly reduces the number of Gribov copies but un- 
fortunately not to a single one. The Gribov ambiguity is 
therefore still present in this model. In the first imple- 
mentations of the GZ model, the predictions for the prop- 
agators were consistent with the scaling solution. More 
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recently a refined version of this model was introduced 
[23 | , that takes into account the appearance of some con- 
densates. The results are then consistent with the mas- 
sive solution of SD equations. 

The results obtained in the SD approach triggered a 
large activity from the lattice community. Since the pre- 
dictions of the SD approaches are mostly for gauge depen- 
dent quantities, it was necessary to implement a gauge 
fixing procedure in the lattice, at odds with usual simula- 
tions. One of the merits of the precursor work of Alkofer 
and von Smekal Q was to discard very singular so- 
lutions found previously for the gluon propagator and 
this was confirmed by the lattice simulations. However, 
a broad consensus on the details of the simulation re- 
sults in d = 3 and 4 was achieved only recently, mainly 
through large lattice simulations |25l429| . It is now well 
established that in these dimensions, the simulations sup- 
port the massive solution. The d = 2 case is different in 
this respect [l^, H^l and is well described by the scal- 
ing solution and the explanation of this difference is still 
unclear (see however [3l|). 

Recently [s^], we proposed an alternative to the an- 
alytical methods described above. Its motivation relies 
on the fact that we do not have so far a fully justified 
covariant gauge fixing procedure that can be handled in 
analytical calculations. Ideally, we should try to con- 
struct such a unambiguous gauge-fixing procedure, but 
this is an extremely hard task. Alternatively, we propose 
to construct a model that uses as guiding principle the 
behavior of correlation functions observed on the lattice 
and that respects as many important known properties 
of correlation functions as possible. More precisely, since 
the gluon propagator is observed to be massive in the lat- 
tice simulations (for d > 2), we propose to include this 
mass term directly at the bare level. (This corresponds 
to a particular case of the Curci- Ferrari model ^33,].) This 
modifies the theory in the IR but preserves the standard 
FP predictions for momenta p ^ m at all orders of per- 
turbation theory. In particular, the mass term does not 
spoil the renormalizability. This model reproduces at one 
loop order with excellent precision the lattice predictions. 
We also showed that the spectral function of the gluons 
is not positive-definite, in agreement with other studies 

Our aims in the present article are the following. First, 
we give a detailed version of the calculations presented 
in [32l |. including some other one loop calculations that 
were not presented in the previous work for lack of space. 
Second, we show that an appropriate renormalization 
scheme can be chosen that does not present a Landau 
pole in the IR. Third, we discuss in detail the size of 
higher loop corrections showing that, in fact, they seem 
to be rather small. This opens the door to perturba- 
tive calculations in QCD in the IR regime, with a huge 
number of possible applications. Fourth, we analyze the 
d — 2 case and explain why in the present model it is 
very different from the d > 2 case. 

The article is organized as follows. In Sect. II, we 



describe in detail the model and its general properties. In 
Sect. Ill, we present the 1-loop perturbative calculation 
of the 2-point functions and compare our results with 
the lattice data in d = 4 and d — 3. In Sect. IV, we 
perform a renormalization group (RG) analysis of these 
propagators. We propose two different renormalization 
schemes, one of which is shown to be IR safe in the sense 
that it leads to no Landau pole in the IR. We compare 
the RG results with the lattice data in d = 4 and d = 3. 
In Sect. V, we discuss the d = 2 case. Some technical 
details are presented in two appendices. 



II. THE MODEL 

As said in the Introduction, the FP action is not jus- 
tified at a non-perturbative level. In this section, we 
present a modification of the FP action in the Landau 
gauge, based on phenomenological considerations. Our 
main guide is the observation that the gluon propaga- 
tor tends to a finite positive value in the IR for d > 2. 
We propose to impose this property at the tree level by 
adding a mass term for the gluon in the FP action. We 
do not change the ghost sector since the ghost propagator 
is found to be IR divergent in the simulations. This may 
also be motivated by assuming that the shift symmetry 
c — > c -f est. is preserved beyond perturbation theory. 
Moreover, we do not want to modify interactions in the 
action so as to preserve the ultraviolet (UV) behavior of 
the theory and maintain the predictions of perturbative 
QCD (or gluodynamics) for momenta much larger than 
Aqc d ■ From this analysis it is clear that if we choose not 
to modify the field content of the theory, mass terms are 
the only local and renormalizable modifications of the FP 
action that do not affect the UV behavior of the model. 
If we restrict to local terms, the only other possible way 
to modify the action is to introduce new fields as done in 
the GZ model [11, [2l|, [22| . The origin of the appearance 
of the effective mass term in the IR is a complex prob- 
lem and could be related with the Gribov ambigui t y, with 
some kind of condensates (see for example [H, [H, [s^l ) , or 
with other nonperturbative effects, but our phenomeno- 
logical approach does not rely on which of these scenarios 
is valid. 

This analysis leads us to consider the Landau-gauge 
FP euclidean Lagrangian for pure gluodynamics, supple- 
mented with a gluon mass term: 

^ = ^iPP' + d,c%D,cr + ih'^d.Al + — (Ap2, (1) 

where {D^,cY = d^^c" + gf^^A^^c", the field strength 
F;^^ = d^,Al - d^Al + gf^^^A^^Al are expressed in terms 
of the coupling constant g and the latin indices corre- 
spond to the SU{N) gauge group. The Lagrangian ([T]) 
corresponds to a particular case of the Curci-Ferrari 
model [33]. At the tree level, the gluon propagator is 
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massive and transverse in momentum space: 

1 



(2) 



with P^i(p) 



p^Pv/v^- It is interesting to note 



that the spectral density associated with the propaga- 
tor ^ is positive and therefore there is no violation of 
positivity at the tree level. As discussed in |32i], viola- 
tions of positivity are present in the model but they are 
caused by fluctuations. 

The gluon propagator observed in the lattice is not 
compatible at a quantitative level with the bare propa- 
gator ^ and we will show below that, by including the 
one-loop corrections, one obtains propagators for gluons 
and ghosts that are in impressive agreement with those 
obtained in the lattice in d = 4 and d = 3. 

Let us mention that a mass term has been used to im- 
prove perturbative QCD results in order to reproduce the 
phenomenology of Strong Interactions |38-l40l . Moreover, 
there are successful confinement models 4l| that use ac- 
tions including a gluon mass term. The difference with 
the model used in those works is that the Curci-Ferrari 
model is renormalizable, allowing to perform perturba- 
tive calculations at any order. This also implies that the 
UV beta function identifies with the standard results of 
the FP procedure. 

When analyzing the model described above, we must 
face the problem that the mass term breaks the BRST 
symmetry (42l . l43j which is very important in the pertur- 
bative analysis. This symmetry has the form: 



bc„b „c 



(3) 



where is a global grassmanian parameter. The BRST 
symmetry is in general used to prove the renormalizabil- 
ity of the theory. However, the breaking of the BRST 
symmetry by the mass term is soft and therefore does 
not spoil renormalizability [s^ |43| . 

The BRST symmetry is also used to reduce the state 
space to the physical space, in which the theory is uni- 
tary (at least at the perturbative level) and breaking this 
symmetry spoils the standard proof of unitarity. This 
problem is actually common to essentially all methods 
that try to go beyond the standard perturbation theory 
(as the GZ model) because they all break the standard 
BRST symmetry. In this respect, the model considered 
here is not in a worse position than other approaches 
considered in the field. We must insist that this model 
is equivalent to the standard FP model in the UV limit 
p ^ Aqc£) if m ^ ^QCD- This means that in the domain 
of validity of standard perturbation theory, the model is 
as unitary as QCD. The unitarity of the model in other 
momentum regimes is of course an important open prob- 
lem, as it is in all gauge fixings in which standard BRST 
symmetry is broken. 

The model with Lagrangian ([l}, as a particular case 
of the Curci-Ferrari model, has a pseudo-BRST symme- 





FIG. 1: First line: three diagrams contributing to the gluon 
self-energy. Second line: diagram contributing to the ghost 
self-energy 



try (not nilpotent) that has the same form as the stan- 
dard BRST ([3]) except for the h variation which reads 
Sih"" = rj rri^c"'. On top of this symmetry, the Lagrangian 
has all the standard symmetries of the FP action for 
the Landau gauge. This includes the shift in antighost 
c — )■ c + cst., a symplectic group [i^, and four gauged su- 
persymmetries [53| HE] • As a consequence, the mass [13] 
and coupling constant (48j renormalization factors are 
fixed in terms of gluon and ghost field renormalizations. 
More precisely, by using the standard definition of renor- 
malization factors (here the subscript 'B' denote bare 
quantities corresponding to their respective renormalized 
quantities without subscript): 



ZaA 



a ^ 



a / '7 a —a / y -a 

= Zgg 



2 ry 1 

nin = Z„,2m 



(4) 



one can prove that the divergent part of the renormal- 
ization factors (or, similarly, the renormalization fac- 
tors themselves in a MS scheme) verify the two non- 
renormalization theorems 143, : 



(5) 
(6) 



ZaZc = 1, 



These non-renormalization theorems are particular cases 
of those found in the Curci-Ferrari model for any gauge 
parameter ^, as proven recently [i^ll^. The finite parts 
of renormalized vertices also verify non-renormalization 
theorems, as will be discussed below, but their explicit 
form depend on the considered renormalization scheme. 



III. PERTURBATIVE CALCULATION OF 
2-POINT CORRELATION FUNCTIONS 

We present in this section a strict 1-loop calculation 
(that is, without taking into account RG effects) of gluon 
and ghost propagators in arbitrary dimension. These cor- 
relators require the calculation of four Feynman diagrams 
as shown in Fig. [TJ The corresponding results in d = 4 
and d = 3 are then compared with lattice results, show- 
ing a very good agreement (for momenta not much larger 
than to). Some of these results were presented in [s^j but 
for completeness we review all of them here. 
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Observing that the gluon propagator is transverse in 
Landau gauge, it is convenient to parametrize the gluon 
G°fl{p) and the ghost G'^''{p) propagators in the form: 



The remaining integral can be performed analytically in 
integer dimensions as discussed below. Let us only men- 
tion here the UV divergent part when d — 7> 4: 



G'^'ip) = 5abF{p)/p\ G^lip) = P^Ap)5abG{p). (7) 

The F{p) is known as the ghost dressing function and the 
scalar function G(p) will be referred to as the gluon prop- 
agator below. We also define the 2-point vertex functions 



T'^i\p)^G-\p), 

Tf;^{p)^p^F-\p). 



(8) 



We choose the following renormalization conditions for 
2-point functions: 



T'^i\p = 0) = m^ T'^i\p ^^i)^m^+^? 



(9) 



We use the Taylor scheme for the coupling [48j], defin- 
ing the coupling constant from the ghost-antighost-gluon 
vertex when the ghost momentum is zero. Using the fact 
that this vertex has no quantum corrections, one deduces 
that ^ is valid also for the finite part of the renormal- 
ization factors. 



r 



(2),1 loop 



g^Np^ 



(11) 



2 167r2(4-d)' 

This is in agreement with previous calculations in the 
literature H-ISOll. 



Let us now consider the 2-point vertex for gluons. The 
first diagram contributing to this vertex function is the 
gluon tadpole (first diagram in the first line of Fig. [1]) 
and gives a contribution: 



.(2) 



d~2 



(47r)'i/2 d 



r(i-x). 



(12) 



A. Strict perturbative results for 2-point 
correlation functions 

Let us first consider the calculation of the ghost self- 
energy in arbitrary dimension using dimensional regular- 
ization. Only one diagram contributes (see second line of 
Fig. [T]). The only difference with respect to standard 
Yang-Mills calculations is the form of the bare gluon 
propagator ([2|). Introducing Feynman parameters and 
performing internal momentum integrals one easily ar- 
rives at the following expression for the ghost 2-point 
vertex: 



r 



(2),1 loop 



X / dx{{xm^ +x{l- x)p^Yf^~^ 

2 xra^ + x{l — x)p'^ ~ 
^ 2-d . 

x{l — x)p^ \ 



X [m^ +p^(l - xf 
~{x{l~x)p^f'^-^{fx^ + 



2-d 



(10) 



The only effect of this diagram is to renormalize the mass. 
The second term contributing to the gluon self-energy is 
the ghost sunset (second diagram in first line of Fig. [T]) 
which reads: 



.(2) 
A, 2 



{P) 



(r(|))V(. 



(47r)d/2 



(2 - d)T{d) 



(13) 



Finally, the third diagram contributing to the gluon self- 
energy is the gluon sunset (third diagram in the first line 
of Fig. [1]). Again, introducing Feynman parameters, we 
obtain: 
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d'^q j 2 (x^p'^d + (a;2 + 2) q^p^ + q^) ( 1 1 



m?'! \ (q2 __|^ + x{l — x)p'^)'^ [rn^ + q^ + x{l — x)p'^)^ 



2p2 + g2 _^ (d + 2)g6 + (8x2 _ + 2d (.t^ - a; + l) + 5) gV + (d + 2) (x^ - a; + l) gV 



i4 



(m2 + 92+;^(l_2;)p2)2 d(d + 2)m^ 

/ 1 1 11 



(m2(l — x) + q2 -|- 2;(1 — x)p'^)^ [mp- + q^ + x{l — x)p^)^ {xmn^ + q^ + x{l — x)p'^)^ [q^ + xil — x)p^)'^ 



(14) 



The integral in the internal momentum q can be done an- 
alytically in arbitrary dimensions but the result is lengthy 
and not particularly illuminating. We can compute the 
divergent part of the gluon 2-point function, which reads: 



-r^(2),l loop/ X d^i 



167r2(4-d) V2 



13 



-P 



(15) 



Together with (llip , this leads to the determination of the 
divergent part of the renormalization factors at one loop: 



35 g^N 1 



Zr 



Za = 1 



6 167r2 4 
3 g^N 1 
2 167r2 4 - d 
13 g^N 1 



d ' 



(16) 



3 167r2 4 - d ' 
which coincide with previous results [13, [13] • 

B. Results in d = 4 

When d — 4 the expressions for correlation functions 
(|10I12I13I14P diverge in the UV. In order to calculate the 
renormalizcd 2-point functions, one must consider the 
sum of those self-energies with a bare contribution with 
renormalization factors. When d approaches an integer 
dimension, the remaining integrals can be performed an- 
alytically. 

Once the renormalization conditions are imposed (see 
Eq. ini the renormalizcd functions F{p) and G{p) are fi- 
nite in the limit d — )■ 4 and read: 



G-\p)/m^ = s + l + 



g^Ns 



nis-^ - 2s-2 



3847r2 . 

+ (2 - s2) log s 2(s-i + I f (s2 - lOs + 1) log(l + s) 

V4 + s-^/i' 



+ (4s-i + 1)3/2 |-^2 _ 20s + 12) log 



\/4 + s + 



F-i(p) = 1 + |^{ - slogs + (s + l)3s-2 log(s + 1) 



(17) 
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FIG. 2: Four-dimensional correlation functions for 5(7(2) 
gauge group. The results of the present work (red curve) are 
compared with lattice data of [2g | (blue points). Top figure: 
gluon propagator. Bottom figure: ghost dressing function. 



where s — p^ /m? . 

In Fig. [21 we compare these expressions for the SU (2) 
gauge group with the lattice simulations of [25|. The 
best choice of parameter is g = 7.5 and m = 0.68 GeV 
when the normalization prescriptions are imposed at 
H — 1 GeV. One observes that both gluon and ghost 
propagators can be fitted with the same choice of param- 
eters in a very satisfactory way. Note that the normal- 
ization conditions used in lattice simulations are not ([9]) . 
Accordingly, it is necessary to introduce a global mul- 
tiplicative renormalization factor when comparing the 
curves. 

We have also compared our results with the data of 
two different lattice studies [1^, [2^ for the SU{?)) group. 
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3 4 5 

p (GeV) 

FIG. 3: Four-dimensional correlation functions for SU{3) 
gauge group. The results of strict perturbation theory (red 
solid curve) are compared with lattice data of [23] (blue open 
circles) and [2^ (blue crosses). The (black dashed) curve 
is obtained by the zero momentum prescription RG (see 
Sect. IIV A|) . Top figure: gluon propagator. Middle figure: 
gluon propagator times . Bottom figure: ghost dressing 
function. 



der to make visible the UV regime. The best choice of 
parameters is g = 4.9 and m = 0.54 GeV (again with 
the renormahzation prescription imposed at = 1 GeV) 
and it leads to a very satisfying agreement for momenta 
p < 2 GeV. Beyond 2 GeV, the agreement is not as good 
but this is not a surprise because expressions (jl7p are 
1-loop results obtained from a fixed coupling constant 
calculation in a fixed renormahzation point. It is well- 
known that in order to analyze the regime p 3> to, one 
must take into account RG effects and in particular the 
running of the coupling. The corresponding procedure is 
presented in detail in the next sections. Once RG effects 
are included, the agreement is essentially within error 
bars for p > m as is also shown in Fig. [3] In any case, it 
is obvious that when p ^ m, the model ([I} reproduces 
correctly the high momentum regime once RG effects are 
taken into account because it just coincides in that case 
with the standard FP model. 

A natural question to raise at this point is whether 
the parameters found by the previous fitting procedure 
are compatible with other determinations. We choose 
to compare g with the results of [5l| which make use of 
a very similar renormahzation scheme: the coupling is 
defined through the Taylor scheme and the ghost renor- 
mahzation factor is fixed by the second line of Eq. 
The only difference is that the gluon renormahzation fac- 
tor is fixed by the condition r^''(p = /x) = /x^. Because 
of the non-renormalization theorem of the coupling, this 
amounts to a relation between the coupling constant g^^^ 
of (SlJ and our which reads: 



(18) 



We extract from Fig. 1 of Ref. the value g^^^ ~ 3.5 at 
^ =1 GeV while the RHS of ^ leads to 4.3. The error 
is about 20% which, as discussed below is the typical 
estimate of higher loop corrections, see Sect IIV El 

An interesting feature of the 1-loop gluon propaga- 
tor is that it is increasing in the IR. In fact, the in- 
verse propagator behaves at small momenta as m? -\- 
7V.gV/(1927r2)log(pVTO2) + 0{p^). This prediction of 
our calculation has a very small effect for d = 4 and it is 
not visible in Figs. I2I3I but appears clearly in d = 3, as 
shown below. 



These two data sets have different overall momentum 
scale and we have rescaled the momenta of the data of [2^ 
for superimposing them with those of [29} . Contrarily to 
the SU (2) case where data are only available in a small 
momentum interval up to 1.9 GcV, the data for S'[/(3) 
were computed for momenta up to 8 GcV. Consequently 
one can explore the cross-over from standard perturba- 
tive results when p is significantly larger than m (requir- 
ing RG methods) and an IR regime similar to the one 
analyzed in the SU{2) case. We represent in Fig. |3] the 
gluon propagator and the ghost dressing function. We 
present also the gluon dressing function p^ G{p) in or- 



C. Results in d = 3 

Let us now consider the three-dimensional case. As in 
d = 4, the 2-point vertex functions can be obtained in an 
explicit form: 



^(2), Hoop. ^ g'^mN f ^ 



2(s-M)^arctan(^/s) 



(19) 
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r:f'-(rt = jg^(-4(5.^ + 7.-i)^ 

+ TT (s^ - 2) - 4(.s + 1)2 (s^ - 6s + 1) arctan (v^) 
+ 2s(s + 4) (s^ _ 12s + 8) arctan ) , 

(20) 

where again s = /m^. 

These expressions are finite because in d = 4 all di- 
vergences are logarithmic and lowering the dimension re- 
duces the level of UV divergence of any diagram. Ac- 
cordingly, a possible renormalization scheme could be to 
use as renormalized parameters just the bare ones (that 
are finite) . However it is well known that is is more con- 
venient to chose a renormalization scheme where renor- 
malized parameters are defined at a running scale. This 
corresponds to a finite renormalization with respect to 
the bare parameters. This is convenient in order to avoid 
large perturbative corrections. Moreover, choosing renor- 
malized parameters at a running scale will allow us to 
perform a RG analysis in the next sections (that improves 
considerably the results). In the present section we use 
the renormalization conditions © • 

The model ^ is able to account for the main features 
of giuon and ghost propagators found in lattice simula- 
tions. In Fig. |4] the results of our calculation with the 
best fit parameters g = 3.7 VGeV and m ~ 0.89 GeV 
for /X = 1 GeV are compared with d = 3 simulations per- 
formed with the gauge group SU (2) [2^ . We observe that 
the best fit for gluon and ghost propagators are not as 
good as in d = 4. We will show below that the inclusion 
of RG effects improve significantly one-loop results. In 
any case, our calculation reproduces the finite IR gluon 
propagator and ghost dressing function. It also repro- 
duces the increasing behavior of the gluon propagator in 
the IR. Indeed, an expansion of the inverse propagator 
at low momentum leads to — Ng^p/QA + 0{p'^). 

IV. RENORMALIZATION GROUP ANALYSIS 

We showed in the previous section that the 1-loop per- 
turbative results compare very well with the lattice data 
up to energies of order 2 GeV and get worse at higher 
energies. This is not a surprise since in d = 4 the per- 
turbation theory generates logarithms of the momentum 
divided by some momentum scale (the mass or the renor- 
malization point /i). The pure perturbation theory fails 
at energies higher than 2 GeV because these logarithms 
are large and one needs to implement the ideas of RG. 
Actually since we have massless modes (ghosts) in the 
theory, it is expected a priori that the RG is necessary 
also for p <C TO. We will see however that in this par- 
ticular model, the IR behavior is milder than expected, 
essentially because ghosts interact by the exchanging glu- 
ons which are massive. 

For d < 4, the UV behavior does not present large 
logarithms, and consequently the use of the RG is not 
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FIG. 4: Three-dimensional correlation functions for SU{2) 
gauge group. The results of strict perturbation theory (red 
solid curve) are compared with lattice data of (blue 
points) . The (black dashed) curve is the IR-safe RG improved 
result (see Sect. IIV B[l . Top figure: gluon propagator. Bot- 
tom figure: ghost dressing function. 



mandatory in that regime. However, in practice, taking 
into account RG effects may improve the quantitative re- 
sults. This idea is natural, but we are not able to test it 
because there is no UV data available in d = 3 or d = 2. 
For d < 4, a RG adapted for the IR is much more im- 
portant than for d = 4 and at the end of the present 
section we show that such a RG procedure improves con- 
siderably the results in d = 3. The comprehension of the 
d — 2 case also improves considerably by the inclusion of 
RG effects but this discussion is postponed to Sect. |Vl 



A. Vanishing-momentum prescription scheme in 

d = 4 



We recall the main steps of renormalization, mainly 
to fix the notations. We concentrate in this section on 
the renormalization prescriptions that were described in 
Sect, ini We define the (3 functions and anomalous di- 
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mensions as: 



d9_ 

dm? 



dyL 
dlogZA 



dfi 
d log Zc 



dfi 



(21) 
(22) 
(23) 
(24) 



We can then use the RG equation: 
to relate the 2-point vertex functions at different scales: 



(25) 



= ZA{^^)TAAiPl^J■o,9i^J■o),m'^{^J■a)), (26) 
r& (p,Ai,5(A*),"i^(M)) 

= ZcifJ-)T'i.\p,fj.o,g{iio),m^{fJ-o))- (27) 

where g{p) and m?{^,) are obtained by integration of the 
beta functions with initial conditions given at some scale 
/io and: 



Jfj.0 A* 

log Zcin) = f ^^7, {g{^i'), m^{n')) . 



(28) 



Using the normalization conditions described in Sect. HIl 
we then get: 



rV(p,Mo,ff,"T. 



Zcip) 



(29) 



(30) 



The renormalization scheme considered here leads to 
a relation between the /3 functions and the anomalous 
dimensions 7a and 7c. As explained before, in the Taylor 
scheme the non-renormalization theorem ([5]) is also valid 
for the finite parts. This leads to the following expression 
for the /? function of the coupling constant: 



f lA{9,m^) 
V 2 



+ 7c(5,™^) 



(31) 



The choice of a normalization prescription at vanishing 
momentum for the gluon 2-point vertex simplifies greatly 



the analytical treatment because as shown in the Ap- 
pendix |X] the following relation is satisfied for the finite 
parts of the propagators and dressing functions: 



vf{p = Q) = ZAmlFB{p^Q). 
It is then easy to prove that: 

/3m2(5,»7i^) =m^-iA{g,m^)- 



(32) 



(33) 



The perturbative results together with the renormal- 
ization scheme (j9l) enables us to compute the anomalous 
dimensions. We get: 



lA 



g'N 



t{3At^ ~ 175t + 6) ~ 2t^ logt 



+ 2{t + l f{2t^ - nt^ + 20t - 3) log(i -I- 1) 
-I- 2t^/'^VTT4: {t^ - 9t^ + 20t - 36) 
+ 4- Vi> 



log( 



g'N 



-5^(2(t+l)t-i3iogi 

+ (< + l)2(t-2)log(i + l)), 



(35) 



where t = fi^/rn^. 

These expressions for the anomalous dimensions lead, 
together with Eqs. ([5T|) and p3l) . to the expressions of the 
f3 functions. Observe that these /3 functions are not only 
function of the coupling constant g (as for example if we 
use the MS scheme), but also depend on the dimension- 
less ratio /.t^/m^. In principle, once these flow equations 
are integrated with some appropriate initial conditions, 
we would need to make another integral to determine 
za(m) and Zc(m) (see eq. (|28l)) in order to determine the 
inverse propagators. However, since the flows of g and 
are expressed in terms of the anomalous dimensions 
for the ghost and gluon flelds, these integrals can be done 
explicitly: 



^a(/^) = 



TO^(/i) 

m2(/zo)' 



5(Mo) V '^^(m) ' 



(36) 



(37) 



Studying the behavior of these anomalous dimensions 
for /i much larger and much smaller than m together with 
eq. (I5T|) . we find: 



Pg 



£N_ll 

W-rr^ 3 
/AT J_ 
T&Tr^ 12 



if 
if 



1-1 » m, 
fi <C m. 



(38) 



In the UV {fi ^ m) we obtain the standard, univer- 
sal, /3 function. This is in contrast to other approaches 
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where masses are introduced in Yang-Mills theory but 



411. In the 



modifying the UV behavior of the model 
IR (/U <C m), we observe that, although smaller in ab- 
solute value than in the UV, /3g remains negative which 
means that the coupling constant diverges at a finite en- 
ergy scale: we encounter a Landau pole in the integra- 
tion. 

However, we stress that there are no large quantum 
corrections when p ^ m (at least when d > 2) because 
there is no IR divergences. In this configuration, a natu- 
ral scheme is to use strict perturbation theory for p m 
and a RG improvement when p ^ m. The absence of 
IR divergences might be surprising since the ghosts are 
massless. A simple way of understanding this is to ob- 
serve that all diagrams, except those with a single ghost 
loop, include at least one gluon propagator, which reg- 
ularize the IR behavior (this issue is discussed in more 
detail in Appendix [B|). In this sense a purely pertur- 
bative calculation for momenta p ^ 1 GeV and a RG 
improvement for UV regime {p ^ 1 GeV) is fully justi- 
fied. In practice, since the RG flow presents a Landau 
pole at an energy scale not far from IGeV, we want to 
make the matching at a higher energy, but lower than 2 
GeV where strict perturbation theory begins to become 
untrustable. A compromise is to perform the matching 
condition at 1.5 GeV and this gives impressively good 
results, as can be seen in Fig. [31 



B. IR safe scheme in d = 4 

A drawback of the procedure described in the previous 
section is that there is some arbitrariness on the choice 
of the matching scale to join strict perturbative and RG 
results. In this sense a unified procedure, valid at all 
scales would be more satisfactory and this motivates the 
search for an IR safe RG scheme. 

In this respect, it is interesting to note that the IR 
Landau pole that appears in the previous RG scheme 
is not a physical effect but is actually a consequence of 
our renormalization scheme Indeed, the prescrip- 

tions we chose there are not consistent, if fi is too small, 
with the 1-loop results (and with the actual behavior ob- 
served in lattice simulations in d = 3) that lead to an 
increasing propagator in the IR (see Sect. IIIll) . In order 
to avoid this problem and be able to take renormaliza- 
tion prescriptions at arbitrary scales, we present here a 
scheme which is not based on an explicit constraint of 
the 2-points vertex function but which comes from the 
nonrenormalization theorem (jH]). In practice, we replace 
the first condition in (|9]) by imposing the relation: 



1, 



(39) 



not only for the divergent parts but also for the finite 
parts. The renormalized mass does not correspond 
any more to the value of the renormalized propagator at 
zero momentum and accordingly there is no contradiction 
with the increase of the propagator in the IR. 



Within this new scheme, few formulas must be modi- 
fied. The /3 function for the mass is now changed to: 

P„,2{g,m^) ^ m'^{-/A{g,m'^) +"fc{g,m'^)). (40) 

The anomalous dimension for the ghost is unchanged (see 
eq. ([55]) ) while the anomalous dimension for the gluon 
field reads: 

+ {-m^ + 7 At - I2)t + t^ log(i) 
- t^l'^^/tTl {t^ - 9t^ + 20t - 36) 



(41) 



X log 



where, as before, t = /i^/m^. Finally, the expressions 
relating Zj^lp) and Zc(/i) to the coupling constants read 
now: 



za{p) 



g^ifi) m^(^o) 
g'^iHo) "i2(^) ■ 



(42) 



(43) 



The UV universal behavior of the /? function (for ^ ^ 
m) is unchanged as expected. However, we observe that 
the beta function for the coupling constant in the IR is 
now positive: 



167r2 6 



if 



(44) 



We need to insist on the fact that the universality of the 
beta function is only valid for mass-independent schemes 
or for /i 3> TO up to two-loop order. In particular, for 
/i < TO there is no reason for the beta function to be 
scheme independent and it is typically not the case. The 
result (|44l) means that in the present scheme, the IR be- 
havior is possibly safe. Actually, depending on the initial 
conditions of the flow we do, or do not, hit a Landau pole. 
(Obviously, if we start with a vanishing mass, we retrieve 
the standard 1-loop calculation and we do hit the Lan- 
dau pole, so a sufficiently large initial mass is necessary 
to be IR safe). 

A second consequence of this IR behavior is that if 
the flow does not hit a Landau pole, then the coupling 
constant is attracted towards zero in the IR: the mas- 
sive gaussian fixed point is attractive. This is a very 
appealing property because this justifies the use of per- 
turbation theory in the regime p <^ m, although its use 
at intermediate momenta (p to) is more delicate. The 
validity of perturbation theory at intermediate momenta 
is discussed in Sect. I IV El 

In Fig. [5] the results for the correlation functions ob- 
tained from the numerical integration of the RG equa- 
tions for the coupling and mass is presented. The best 
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FIG. 5: Four-dimensional correlation functions for SU{3) 
gauge group. The results obtained by the IR safe RG (see 
Subsec. IIVB|I (red solid line) are compared with lattice data 
of [i^ (blue open circles) and _2g] (blue crosses). Top figure: 
gluon propagator. Middle figure: gluon propagator times . 
Bottom figure: ghost dressing function. 



choice of parameters at the scale n — 1 GeV is 5 = 3.7 and 
m — 0.39 GeV. One observes a very good agreement of 
the 1-loop RG result when compared with lattice simula- 
tions, particularly, as expected, in the UV and in the IR. 
The only region where a significative departure is seen is 
for intermediate momenta (1 GeV < p < 2 GeV) for the 
gluon dressing function. The results, however, are not as 
good as those obtained with the scheme of Sect. IIV Al 
This may be surprising and we now explain why it is 
not so. First, as already mentioned and explained in 
detail in the Appendix [B] there are no IR divergences 
and consequently the use of RG is not mandatory for 
p < m. Second, the results of the previous section have 
been obtained by a mixed treatment, by using the best 



fit of strict perturbation theory and imposing a match- 
ing with RG results in the UV. The results discussed 
here are obtained with a unified treatment, valid at all 
scales. This is a more challenging situation because we 
do not have the freedom of choosing the matching scale. 
In this sense the IR safe scheme gives a more satisfying 
construction. Third, as discussed in Sect. IIV El even if 
the relevant coupling in the present model is never large, 
it is nevertheless not very small. Accordingly the result 
obtained in the present section are quite satisfactory and 
the surprise is more the level of precision achieved in the 
previous one. A possible explanation is that the choice 
of RG prescriptions that can be directly read from cor- 
relation functions (and not from an implicit RG scheme 
as done in (j39p ). may in practice reduce the contribu- 
tions of higher loop corrections. In any case we may 
expect that two-loop corrections should improve consid- 
erably the present results and correspondingly reduce the 
scheme dependence of them. Finally, the coupling con- 
stant compares much better with the one extracted from 
lattice simulations (see Eq. (|18l) '). Taking into acount 
the difference of renormalization scheme, we now find, at 
/i = 1 GeV, g^'-* = 3.4 which compares very well with the 
lattice data 3.5 [Slf. see Eg. pS]) . 

C. The IR behavior in arbitrary dimensions 

Before discussing the results in d = 3, it is convenient 
to study the IR {fj, ^ in) behavior of RG functions of 
the model in arbitrary dimensions. We saw in the pre- 
vious subsection, that in that regime, the beta function 
for the coupling constant is positive and accordingly the 
gaussian IR fixed point is attractive. When d 7^ 4, the 
coupling constant g has dimension [energy] ^^"'^^/^ and it 
is convenient to introduce the dimensionless coupling 

= (45) 

in order to analyze a process with characteristic momen- 
tum fj,. Of course, it may be appropriate to use powers 
of the mass in order to define a dimensionless coupling. 
This point is discussed below (see Sect. IIVE[) . but for 
the moment let us concentrate on the choice pS)) which 
is the natural definition of the dimensionless coupling in 
the UV regime m). The /3- function associated with 

9 is: 

For dimensional reasons, the function /3g is only a func- 
tion of g and of the dimensionless squared mass: 

^ m^fi^^. (47) 

For any d < 4, as well known, in the UV regime {fi 3> m 
or, equivalently ^ 1), the coupling g approaches 
rapidly a gaussian fixed point because of the dominance 
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of the dimensional term — ^^g. The case li = 4 is also 
dominated by a gaussian fixed point in the UV, because 
of the negative sign of the ff^ term. At intermediate mo- 
menta /i ~ m, the expressions of the /3g function is very 
complicated in arbitrary dimensions and its (non illumi- 
nating) expression involves special functions. However in 
the IR regime fi m things become easier again, and 
one can give an explicit expression for the (3g function at 
one-loop: 



-9 



(4-d)(d-2)7r' 



2 " ■ 22rf+ir (i^±i) sin (^^) 



N~g'. (48) 



The coefficient of the term is regular (and non zero) 
in all dimensions d > 0. One observes the appearence of 
an IR fixed point g: 



N~g.. 



4^ 



:^)sin( 



(4-rf)7r 
2 



(d- 2)7r' 



(49) 



When d — 4, the expression (|48)) approaches the d = 
4 result (|44)) and the fixed point value approaches the 
gaussian limit = 0. 

The coherence of these results relies on the hypothesis 
that fn} — >■ Qo when /i — >■ 0. To check this hypothesis, we 
write the /3 function for the dimensionless mass rr?\ 



d(rn} [1 



2m^ + fj.-^(3^2. (50) 



As before, the UV regime (where g — 0) is compatible 



with ; 



and this can already be seen in the ^ —2m 



term. The intermediate momentum regime (/i ~ to) gives 
complicated expressions, but the IR regime again simpli- 
fies to: 



-2tf+iv^w .;^-y:^!c:^ , (51) 



4JT(!!±i)sia(ii^) 



Note that the pre-factor of g^rri^ of this expression is 
twice the pre-factor of g^ in the eq. (US]). Consequently 
when g^ approaches its IR fixed point: 



Pfn^ ^ (2 - d)rh 
For d > 2, this implies that a regime with jfi^ 



(52) 
DO and 

g ^ g* exists (under the hypothesis that the flow does 
not hit a Landau pole before approaching it). Note that 
the dimensionless square mass diverges as /i^"'', but the 
dimensionful mass goes to zero for d < 4. This is not in- 
compatible with the fact the gluon propagator is finite in 
the zero momentum limit, because in the present scheme, 
the mass is not defined through the renormalization con- 
dition (|9]). It is easy to show that: 



lim rj^-* [p) — 



glm^iHo) 



"^L5^(Mo)' 



(53) 



where rh [p) ~ rhr^sP ■ The d = 2 case is different. In 
that dimension, the coupling constant would approach a 



fixed point in an hypothetical IR regime, but the lead- 
ing contribution to the flow of to^ is zero and one must 
analyze sub-leading contributions. This important par- 
ticular case is analyzed in Sect. [Vj 



D. IR safe scheme in d = 3 

Having considered the general IR behavior, we now 
repeat the calculation of the Sect. IIVBI in the d — 3 
case and compare the results with the lattice simulations. 
The calculation now leads to the following anomalous 
dimensions: 



■ (^nt^ ( - + 2) + 4\/t (7t^ - 29t -f 15) 



g'N 

- 2t (it^ - + m ~ 96) arctan (Vt/2) 
+ 4 (St" - 4^3 + - 15) arctan (Vi) 



+ 2(t + l)(t - 3) arctan (%/i) 



(54) 



(55) 



where again t = ^ jm^ . We then deduce the /3 functions 
for g and rr? via the non-renormalization theorems pip 
and (PD|) . For some initial conditions the flow shows a 
Landau pole, but for the particular values of couplings 
that do match with numerical simulations, there is no 
Landau pole. The dimensionless running mass TO(/i) in- 
creases when /i decreases and we get an IR regime when 
fi <^ m. As discussed in the previous section, this regime 
is characterized by the ^ <^ m limit of the function /3g 
that for d = 3 is 



/3« 



;;3 



Ng_ 
128 ' 



and the corresponding IR fixed point is 
Ngl = 64. 



(56) 



(57) 



In the Sect. IIV El it is discussed why, in fact, these cou- 
plings are not as large as it seems here. We want here 
only to discuss the corresponding consequences for gluon 
and ghosts propagators. The RG equations can be inte- 
grated numerically and from them, the propagators can 
be obtained as discussed in the d = 4. The correspond- 
ing results are presented in Fig. 01 We note that the 
inclusion of these RG running considerably improve the 
results when compared with those of strict one-loop per- 
turbative results and now compare very well with lattice 
simulations. This is not surprising, because even if the 
present model does not have IR divergences it is very 
usual that RG effects play an important role in field the- 
ories below four dimensions even at moderate values of 
the couplings. 
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E. An estimate of higher loop corrections 

The coupling constants that appear along this work 
are frequently relatively large at intermediate scales /x ~ 
1 GeV. This situation is much better than having a Lan- 
dau pole, but it is anyway problematic for a perturba- 
tive analysis. In this section we argue that although the 
coupling constant is not small, the perturbation theory 
seems to be under control anyway. 

First of all, as is well known |52'], the parameter expan- 
sion is not but, because of angular factors and group 
factors: 



(i-4 



u{p) 



(4^)d/2r(rf/2)- 



(58) 



where p is the momentum scale used to define the cou- 
pling constant. In d ^ 4, g has dimension and 
the power of momenta in the previous formula is neces- 
sary in order to have a dimensionless expansion param- 
eter. The angular factors considerably reduces the size 
of the expansion parameter. For example, for c? = 4 and 
TV = 3, this expansion parameter has a maximum at val- 
ues of order one as shown in Fig. [HI 

There is actually another effect that suppresses the ra- 
diative corrections in the IR as we discuss now in de- 
tail. Note that when momenta are much smaller than the 
gluon mass, all diagrams that include internal gluon lines 
are suppressed by inverse powers of the gluon mass. This 
has two consequences. First, as shown in Appendix [B] 
the present model does not present IR divergences for 
d > 2 at non-exceptional momenta. This is surprising 
at first sight since there are massless modes (ghosts), 
but it can be understood because their interactions are 
mediated by massive gluons. A second consequence is 
that in that regime, the leading contributions are those 
with a minimum number of gluon propagators and we 
can consider an effective theory where the only dynami- 
cal degrees of freedom are the (massless) ghosts. In this 
effective theory, gluons appear only as external sources 
coupled to the ghosts via the bare ghost-gluon vertex 
while the ghosts interact via an effective four-point ver- 
tex that behaves as pimg^ / , where pi and p2 are the 
anti-ghost momenta [5^. This implies that the effective 
expansion parameter is suppressed by p^ jm} where p is 
a typical external momentum. A naive interpolation be- 
tween the UV parameter expansion u given in (|58l) and 
the one relevant for the IR would be: 



g'^(p)Np' 



d-2 



p'^ + m? {4nY/^T{d/2){m^ +p^)' 



(59) 



As shown in Fig. [6] this parameter is indeed rela- 
tively small and this makes a perturbative analysis well 
founded. Actually, in a small region around 1 GeV, the 
expansion parameter is of the order of 0.4. This means 
that 2-loop corrections in this region will give contribu- 
tions of order 0.4^ 15%), which fits with the discrep- 
ancy between the 1-loop results and the lattice results 




fj, (GeV) 

FIG. 6: Four dimensional expansion parameter for 5(7(3) 
(color online). Plain curve (red), u{fi) is the dimensionless 
parameter expansion, valid in the UV. Dashed curve (blue), 
u{ii)ii^ /{fi^ + m^{n)) takes into account the IR freeze-out. 
The parameters are fixed at 1 GeV as in Sect. IIV"B] 
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FIG. 7: Three-dimensional expansion parameter for SU{2) 
(color online). Plain curve (red), u{fi) is the dimensionless 
parameter expansion, valid in the UV. Dashed curve (blue), 
u{ii)ii^ /{fi^ + m^lfi)) takes into account the IR freeze-out. 
The parameters are fixed as in Sect. IIV Dl 



shown in Fig. [S] A 2-loop calculation would give a more 
precise estimate of higher order corrections. 

Similar results apply for other values of N and d. As 
an example, the same curve is presented for N = 2 and 
d = 3 in Fig [T) It is shown that the actual values in the 
IR are slightly larger but remain relatively small. 



V. DISCUSSION OF THE d = 2 CASE 

The explicit expression for propagators can also be ob- 
tained in d — 2. As before, there are no UV divergences 
for the same reasons as for d = 3. The corresponding 
result is: 



t(2),11oop 



(P) 



g^ 
Btt 



(slog(s)-(s + l)log(s + l)), (60) 
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p(2),lloop(p^ 



IOtts 



(-2(s-l)2(s + l)log(s + l) 



\/s(s + 4) log 



Vs + 4 - 



Vs + 4 + 

+ 5(s2-2) log(s)), 



(. - 2)2 

(61) 



where, again s = /m? . 

Within the model, the difference between d = 2 and 
d > 2 which is observed in the lattice (see the introduc- 
tion) appears natural. In d = 2, we find that the gluon 
propagator and ghost dressing function develop logarith- 
mic divergences when p — )■ 0: 



F-\p) 



An 
Anm? 



logm, 
Vm/ 

log(^). 



(62) 



Such divergences exclude the possibility of controlling as 
strict one- loop calculation as was used above for d > 2. 
A proper treatment of the c? = 2 case requires a RG 
approach adapted to the IR regime. 

As shown in the appendix [B] the IR convergence of the 
loop integrals in d = 2 does not improve when increasing 
the order of perturbation theory, as it does for d > 2. 
This is the typical situation at the upper critical dimen- 
sion of a model and it is well known that a RG treatment 
is then necessary to study the IR behavior. We therefore 
apply the IR safe scheme described in Sect. II V Bl in the 
d = 2 case. 

The anomalous dimensions read: 



7A 



SttijH ( - log* + 2(*' - 4) log(l + t) 
+ y^(t-2)(t2 + 4t+12)log(^-^ 



'VtTA + Vi 



Anil' 



log(l + 1) 



(63) 



(64) 



with t — /i^/m^. As before, we can extract from them 
the P functions for the mass and coupling constant. As 
discussed at the end of Sect. lIV C| the leading behavior of 
the flow of the dimensionless mass is vanishing in d = 2. 
In order to study the IR behavior, it is convenient to 
consider the following dimensionless combination: 



m 

whose /3 function takes the form: 



(65) 



NX^ log(l + 1) 

' 8tt 



Observe that this /3 function is always negative and there- 
fore leads to a IR Landau pole. The only way round 
would be that rfi^ = 1/^ would run to zero sufficiently 
fast when fj, decreases. However the study of the flow of 
in the IR is at odds with this hypothesis. Indeed, the 
(3 function for m'^ behaves at small rfi^ as: 



-m^ [2 



NX^ log(t) 



(66) 



In conclusion, the case d = 2 is very different from 
higher dimensions: There is no fixed point and the flow 
runs to a Landau pole. We relate this first property with 
the observation that, in d = 2 and for fj, <C m^, the f3 
functions for g and are proportional (see Eqs. 
and ((5T|) ): 



2 ^9 _ Pm- 



(67) 



t 



A similar situation appears in the random mass Ising 
model |[53| which is characterized by two coupling 
constants. The associated /3 functions are found to be 
proportionnal (in the sense of the previous equation) at 
leading order in e expansion. The two-loop calculation 
[s^ . Issj lifts this degeneracy and yields a fixed point with 
coupling constants ^ ^/e. It would be interesting to study 
if the two-loop contributions lift in the same manner the 
degeneracy of the /3 functions for g and w? . 



VI. CONCLUSION 

We have shown in detail that a particular case of 
the Curci-Ferrari model, motivated by phenomenologi- 
cal considerations reproduces quantitatively several non- 
trivial correlation functions of pure gluodynamics in all 
ranges of momenta. The result is obtained from a very 
simple one-loop perturbative calculation and the surpris- 
ing agreement is justified by estimating the size of higher 
loop corrections. Moreover, we have proposed an IR-safe 
RG scheme that does not show a Landau pole. There- 
fore we have at hands a phenomenological model that 
correctly describes Yang-Mills correlators. These results 
question the common idea that the IR properties of QCD 
are beyond the scope of perturbation theory. 

At the conceptual level, it would be more satisfying to 
be able to calculate the mass parameter from first prin- 
ciples and this is probably a nqn-perturbative and very 
difficult issue (see however [s^, [l^)- However, once the 
existence of this mass is accepted, the rest of the analysis 
seems to be purely perturbative. The determination of 
the mass from first principles would be also interesting 
since it would fix the results in terms of the single pa- 
rameter appearing in pure YM. Another important open 
issue is that of the unitarity of the model. A discussed 
in Sect. II, usual proofs of unitarity do not work in this 
model because the nilponcy of the BRST symmetry is 
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broken. The situation is actually common to all schemes 
beyond the FP procedure, such as the GZ model, but it 
does not unavoidably mean that the model breaks uni- 
tarity. It could be that a procedure permits to reduce 
the state space down to a physical Hilbert space (where 
unitarity is recovered) that would not rely on the nilpo- 
tency of the BRST operator. Devising such a procedure 
is an open issue. 

Beyond these fundamental problems, the evidences 
that the model considered here is a good and simple phe- 
nomenological one opens the way to many applications. 
Let us mention here some of them. 3-point correlation 
functions have been measured on the lattice and their 
calculation is a direct extension of the work presented 
here and would give another test of the ability of the 
model to reproduce YM correlation functions. It would 
also be interesting to compute the 2-loop contributions 
to the ghost and giuon propagators, that would in par- 
ticular give a direct measure of the size of higher order 
corrections. The calculation of RG flows and 2-point cor- 
relation functions can be easily extended at finite tem- 
perature and many thermodynamic properties can be ex- 
tracted from it. Of great physical interest would be to 
calculate within this model the quark-antiquark poten- 
tial in the quenched approximation, although it is hard 
to believe that confinement would appear in such a simple 
analysis. Let us mention finally that introducing matter 
fields is straightforward and studying the influence of the 
quarks on the ghost and gluon correlation functions is a 
natural extension of this work. 
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Appendix A: Non-renormalization theorem for the 
mass 

We derive in this section the non-renormalization the- 
orem for the mass. We used this result in Sect. IIV Al to 
simplify the RG analysis. The main part of the proof 
is done with bare quantities but in order to simplify the 
notations we omit the subscript B. 

The proof makes use of the Slavnov- Taylor identity and 
we therefore introduce sources not only for the primary 
flelds Afj,, c, c and h, but also for their BRST variations. 
The partition function then reads: 



where 



J (fxiJ^A, 



(Al) 



(A2) 



■ cx + + KsA^ + Lsc) 



and S is the integral of the lagrangian ([T]). The BRST 
variations s of the fields are defined as the prefactors of 



rj in the right hand sides of ([3]). It is not necessary to 
introduce sources for the variation of h or for the vari- 
ations of variations of the primary fields since these are 
either vanishing or linear in the primary fields. We also 
introduce the Legendre transform as: 

T + W = j d'^x {Jf,A^ + + Rh) . (A3) 

We first make the change of variables c{x) — >■ c{x) + 
fj{x). This yelds the equation: 



6T 



sr 



6K?^{x) Sc'^ix) 



(A4) 



Deriving this equation with respect to c^{y) and tak- 
ing the Fourier transform (our convention is f{p) = 
/ d'^xexp{ipx)f{x)), we get: 

rfUp) = ~^P,r'^ilJp). (A5) 

We moreover need the Slavnov- Taylor equation, which 
reads: 



d'^3 



sr ST 



sr ST 



SAf^ {x) SKf^ {x) (Sc^ {x) SL°- (x) 



ST 



ST 



(A6) 



= 0. 



We derive this expression with respect to ^(^(y) and c'^{z) 
and take this expression at vanishing sources. Using the 
fact that the vacuum does not break the ghost number 
conservation, only two terms contribute and making the 
Fourier transform we get: 



2„ rfcc 



(A7) 



where we used the fact that the h sector is not renormal- 
ized. 

{'2') 

Now, using the fact that r^„^6(p) is analytic at low 

momentum and behaves as (5"''(5pi,G'^^(0), and contract- 
ing Eq. (|A7p with p^, we obtain at low energy: 

GB\mB\^) = ^l, (A8) 

where we have introduced the subscript B to recall that 
all this calculation was done (as emphasized at the be- 
ginning of the appendix) with bare quantities. 



Appendix B: Infrared behavior of the Feynman 
diagrams 

In this appendix wc study the IR corrections coming 
from higher order diagrams. Our aim is twofold. First we 
check that the regular IR behavior found at 1 loop is not 
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modified by higher order corrections. Second we study 
the IR behavior of the low energy effective theory for the 
ghosts obtained by considering only a minimal number 
of gluon propagators (see Sect. IIV Ep . 

Consider first a generic diagram of the model ([Ij. It 
has Na external gluon legs, Nc external ghost/antighost 
legs, W3 three-gluon vertices, W4 four-gluon vertices and Vc 
ghost-antighost-gluon vertices. Looking at the bare ver- 
tex, we see that each external anti-ghost leg comes with 
an external momentum. In the following when we speak 
of diagrams, it is understood that these trivial external 
momenta dependence have been factorized. 

Let us calculate the superficial degree of IR divergence 
cj of this diagram for vanishing external momenta, which 
characterize the IR behavior of the associated integrals. 
A negative superficial degree of IR divergence indicates 
that the corresponding diagram is IR divergent (at zero 
momenta). Each ghost propagator contributes -2 but the 
gluon propagators being massive do not contribute. The 
3-point interactions being proportional to a momentum 
contribute -1-1 and each loop contributes d. Using text- 
book techniques [11], it is easy to prove that: 

, d d-l d , d-2 , , 
u;>d--NA —N, + -V3 + dvi + ^-v,. (Bl) 

Observing that the pre- factors of W3, W4 and Vc are posi- 
tive in d > 2, we conclude that at fixed number of exter- 
nal legs, increasing the number of vertices suppress the 
contributions in the IR. Moreover we check that the two- 
point vertex functions at zero momenta are finite at all 
loops. The case d < 2 is different because increasing the 
number of ghost-gluon vertices does not improve (or even 
worsen for d < 2) the IR behavior. In this sense, d = 2 



plays the role of an upper critical dimension. 

The situation is, as expected, more favorable if we con- 
sider non-exceptional momentum configurations. Using 
the method described in [5I], we find that the degree of 
IR divergence is: 

d d-l d + 2 , d-2 ,^ , 

w > ^Na + -^Nc + -^V3 + dv4 + Wc, (B2) 

where now Na and Nc are the external soft legs of the 
diagram and U3, V4 and Vc are the number of vertices 
attached only to soft momenta. This expression is always 
positive for d > 2 but again d — 2 plays the role of an 
upper critical dimension. 

We now come to the discussion of the IR behavior of 
the diagrams in the effective theory. We consider first 
diagrams with vanishing external momenta (as before, we 
first factorize the external momenta of the anti- ghosts). 
As discussed in Sect. IIVE| the ghosts interact via a 4- 
point vertex which behaves as the product of momenta of 
the anti-ghosts. Denoting v^^c the number of such vertices 
in a diagram, we find the IR degree of divergence: 

ijj^d-NA- ^^^c + {d- 2)vic. (B3) 
Again, we see that, at fixed number of external legs, in- 
creasing the number of vertices suppress the contribu- 
tions in the IR and check that the two-point vertex func- 
tions at zero momenta remain finite at higher loops. Us- 
ing non-exceptional momentum configurations, we get: 

u = + {d- 2)vic, (B4) 

which ensures that the diagrams are IR finite for d > 2. 
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